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Abstract 

First results of the study of the process e^e~ An by the CMD-2 
collaboration at VEPP-2M are presented for the energy range 1.05-1.38 
GeV. Using an integrated luminosity of 5.8 pb^^ , energy dependence 
of the processes e"'"e~ — > 7r"'"7r~27r'' and e"'"e~ — > 2'k~^2'k~ has been 
measured. Analysis of the diflferential distributions demonstrates the 
dominance of the oitt and ujn intermediate states. Upper limits for the 
contributions of other alternative mechanisms are also placed. 



1 Introduction. 



Investigation of e+e~ annihilation into hadrons at low energies provides 
unique information about interaction of light quarks and spectroscopy of 
their bound states. At present the energy range below J/ijj can not be sat- 
isfactorily described by QCD. Future progress in our understanding of the 
phenomena in this energy range is impossible without accumulation of ex- 
perimental data vitally important to check the predictions of existing theo- 
retical models. In addition, the total cross section of e+e^ annihilation into 
hadrons at low energies as well as the cross sections of exclusive channels 
are necessary for precise calculations of various effects. These include strong 
interaction contributions to vacuum polarization for (g-2)^ and a(M§) Q], 
tests of standard model by the hypothesis of conserved vector current (CVC) 
relating e^e~ — > hadrons to hadronic r-lepton decays ^, determination 
of the QCD parameters based on QCD sum rules Q etc. 

The energy behavior of the total cross section as well as that of the cross 
sections for exclusive channels is complicated and characteristic of various 
broad overlapping resonances (e.g. p, oj, (f> recurrencies) with numerous com- 
mon decay channels having energy thresholds just in this energy range. Pres- 
ence of the broad resonances in the intermediate state necessitates consid- 
eration of the quasistationary states and makes effects of their interference 
important. 

Until recently the investigation of e+e^ annihilation into hadrons was 
restricted by measurements of the cross sections only. Appearance of the 
new detectors with a large solid angle operating at high luminosity colliders 
and providing very large data samples opens qualitatively new possibilities 
for the investigation of the multihadronic production in e~^e~ annihilation. 

Production of four pions is one of the dominant processes of e+e^ an- 
nihilation into hadrons in the energy range from 1.05 to 2.5 GeV. For the 
first time it was observed in Frascati |^ and Novosibirsk First experi- 
ments with limited data samples allowed one to qualitatively study the new 
phenomenon of multiple production of hadrons and estimate the magnitude 
of the corresponding cross sections. Subsequent measurements by different 
groups in Frascati, Orsay and Novosibirsk (see the references in [|l|) provided 
more detailed information on the energy dependence of the cross sections of 
the processes e~^e~ — > 27r"'"27r~ and e^e~ 7r"'"7r~27r'' in comparison with 
the previous measurements. 

One of the main difficulties in the experimental studies of four pion pro- 
duction was caused by the existence of different intermediate states via which 
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the final state could be produced, such as 



e^e OJTT (1) 

e+e" pa (2) 

e+e" ai(1260)7r (3) 

e+e- /ii(1170)7r (4) 

e+e^ p~ (5) 

e+e- a2(1320)7r (6) 

e+e" ^ 7r(1300)7r (7) 



The relative contributions of the above mentioned processes can't be ob- 
tained without the detailed analysis of the process dynamics. First attempts 
of this type were performed by MEA Hand DM1 ||] in the energy range 
above 1.4 GeV and OLYA g and CMD |10) below 1.4 GeV and it was shown 
that the 27r+27r- final state is dominated by the p^vr+Tr" mechanism. ND p4| 
measured the cross section of cjtt production from 1.0 to 1.4 GeV with a mag- 
nitude which was confirmed by the subsequent r decay studies at ARGUS 
as well as by more recent results from CLEG and ALEPH |Q. 

Later the DM2 group tried to perform partial wave analysis (PWA) of the 
mode with four charged pions jl^ in the energy range 1.35 to 2.40 GeV. Their 
analysis was based on the momentum distributions only while the angular 
dependence as well as interference between different waves were not taken into 
account. Although they obtained some evidence for the presence of oi (1260)7r 
and pa states, a mechanism for a substantial part of the cross section was 
not determined. PWA for the mode 7r+7r-27r*' was not performed because of 
the insufficient number of completely reconstructed events. 

The abundance of various possible mechanisms and their complicated in- 
terference results in the necessity of simultaneous analysis of two possible 
final states {2-k^2tt~ and 7r+7r~27r") which requires a general purpose detec- 
tor capable of measuring energies and angles of both charged and neutral 
particles. The first detector of this type operating in the energy range below 
1.4 GeV is the CMD-2 detector at VEPP-2M collider in Novosibirsk |l|]. 
In this work, we present results from a model-dependent analysis of both 
possible channels in e+e" annihilation into four pions based on data col- 
lected with the CMD-2 detector. To describe four pion production we used a 
simple model assuming quasitwoparticle intermediate states and taking into 
account the important effects of the identity of the final pions as well as the 
interference of all possible amplitudes. 
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2 Data sample and event selection 



The analysis described here is based on 5.8 pb~^ of e+e^ data collected at 
center-of-mass energies 2Ei,eam from 1.05 up to 1.38 GeV. The data were 
recorded at the VEPP-2M e+e~ collider of the Budker Institute of Nuclear 
Physics in Novosibirsk, Russia with the CMD-2 detector in 1997. The energy 
range mentioned above was scanned twice with a step of 20 MeV: first by 
increasing energy from 1.05 to 1.37 GeV and then by decreasing energy from 
1.38 to 1.06 GeV. 

The CMD-2 is a general purpose detector consisting of a drift chamber 
(DC) with about 250 /i resolution transverse to the beam and proportional Z- 
chamber used for trigger, both inside a thin (0.4 Xq) superconducting solenoid 
with a field of 1 T. Photons are detected in the barrel Csl calorimeter with 
8-10 % energy resolution and the endcap BGO calorimeter with 6 % energy 
resolution. More details on the detector can be found elsewhere p^ . 

2.1 Event selection 

The present analysis is based on completely reconstructed 7r+7r^27r° events. 
At the initial stage, events with one primary vertex with two opposite sign 
tracks and four or more reconstructed photons were selected. Both tracks 
should come from the interaction region: a distance from the track trajec- 
tories to the beam axis should be less than 0.3 cm and the vertex position 
along the beam axis should be inside ± 10 cm. To reject background from 
collinear events, the acoUinearity angle between tracks in the i? — plane 
should be greater than 0.1 radians. Both tracks are required in the DC fidu- 
cial volume: a 9 angle should be inside 0.54 tt — 0.54 radians. Clusters 
of energy deposition that are not matched with charged track projection are 
paired to form tt^ candidates. These showers must have energies greater than 
20 MeV, and invariant mass of the photon pair must lie within 3a of the tt^ 
mass where a varies between (5 10) MeV. After that a kinematic fit was 
performed assuming the 7r+7r~27r*' hypothesis for all possible tt'^ pairs. For 
further analysis the combination with min(x^) was selected under the con- 
dition /ndf < 2.5. After such selection 22128 events remained in the 
energy range under study. 

2.2 Final event sample 

To understand the dynamics of the process we studied the distribution over 
the recoil mass for one of the tt^'s. Figure |^ shows this distribution at the 
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beam energy of 690 MeV. Each event gives two entries to the histogram 
corresponding to two 7r°. A signal from the ujtt^ final state is clearly seen. 
Points with errors are the data while the smooth line is om' fit corresponding 
to the sum of a Breit-Wigner lo signal convoluted with detector resolution and 
a smooth combinatorial background. The detector contribution to the signal 
width is about 10 MeV. The number of events under the lo peak accounts 
for only ~ 60% of the observed events that indicates at the existence of 
additional intermediate states. 




For further analysis the data sample was subdivided into two classes: 

I. Tmn{\Mrecoa{TT°) - AL\) < 70 MeV 

11. min{\Mrecoui7T'') - MJ) > 70 MeV 

where M^^ is the ui mass. The first class contains mostly W7r° events while their 
admixture in the second class is relatively small, about (1 5)% depending 
on the beam energy. 
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Figs. H and ^ show distributions over Mi„^(7r*7r°) and Mi„t,(7r+7r^) for 
the events in the second class. In the spectrum of Mi„„(7r^7r'^) one can see 
a clear signal of while that for Afi„i,(7r+7r^) is relatively smooth and no 
signal from the p° is observed. The solid lines show our fit including smooth 
combinatorial background and gaussian p signals. Presence of the p^ signal 
and absence of the p° signal lead us to the natural assumption that p mesons 
originate from an isospin 1 resonance. For the case of the resonance with 
7 = one expects production of both neutral and charged p's. 





SOO WOO 



Figure 2: Distribution over 

Minv{'K'^TT^) for 7r+7r~27r° events 
from the class (O) 



Figure 3: Distribution over 

Afim,(7r+7r^) for 7r+7r^27r° events 
from the class (||) 



Possible candidates from the list of the processes (|])-(0) are ai(1260), 
02(1320) and 7r(1300). These resonances have different spin and parity re- 
sulting in different angular distributions for a recoil pion (Table pj). 



Table 1: List of / = 1 resonances with expected angular distributions 









da/dcos{6) 


X^/ndf 


1 


ai(1260) 


1++ 


const 


7.1/7 


2 


7r(1300) 


0-+ 


sin2(6') 


47.3/7 


3 


02(1320) 


2++ 


1 + cos2(6l) 


37.0/7 
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To measure the angular distribution of the p*7r° system (or recoiled tt^) 
we fit the signal in eight ranges of the recoil tt^ angle with respect to 
the beam axis. The measured angular distribution is shown in Fig. ^. Three 
smooth curves show the fits corresponding to mentioned above hypotheses. 
The fit goodness is presented in the last column of Table |^. One can see good 
agreement with the hypothesis of the ai(1260)7r intermediate state. 




Figure 4: Angular distribution of the recoil tt . Solid curves correspond to 
the following hypotheses: 1 - ai(1260)7r; 2 - 7r(1300)7r; 3 - a2(1320)7r 

Thus one can assume that our data can be explained by oj-k^ and ai(1260)7r 
intermediate states although small admixture of other states is not excluded. 



2.3 Fitting method 

For detailed comparison we perform the simulation of the processes (|l|)-(0) 
taking into account their interference as well as the interference of the dia- 
grams differing by permutations of identical final pions (see Appendix ^). 

In order to extract a relative fraction of ai(1260)7r and any other in- 
termediate states, the unbinned maximum likelihood fit has been used p^ . 



Kinematics of each event is completely described by a set of measured four- 
momenta = (e, p) of the final pions with the exception of smearing due 
to detector resolution and radiative effects. This set will be referred to as a 
vector of state Pi = (p7r+iPir-iP7r"jP7r")i where a subscript i stands for the 
event number. 

The theoretical probability density function for Pi depends on process 
dynamics and can be expressed via the matrix element Here 
A stands for a set of unknown model parameters like relative strength and 
phase of interference between the intermediate states. To obtain the optimal 
values of A we minimize fl^ the logarithmic likelihood function 

N2 

LiA)^-Y^lnw{P,,A), (8) 

i=l 

where w{Pi,A) is the probability to observe an event in state Pi, and N2 
is the number of events in the second class. Since the detector resolution 
for the invariant mass of three pions is about 10 MeV, i.e. comparable to 
the u! width, smearing effects are significant for the events of the first class. 
Therefore we use the events of the second class only. To fix the fraction of 
events in the first class we add the following term to the function (||) 

ir{A)-n)^ 

where ri = j^^j^^ is a measured fraction, iVi is the number of events in the 

first class, and r{A) is the expectation of ri. 

The normalized probability density function is expressed via the matrix 
element 

miPuA)\^ 1 



wiP„A) - 



where j is the particle number, and ayis{A) is a model dependent visible 
cross section: 



a,,,{A) = J\m{P,A)f-f[ 



Pi 



The function (Jms{A) is calculated using Monte Carlo technique. For this 
purpose we use a set of events sampled according to the relativistic phase 
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space distribution. In this (A) is calculated as 

-, Nmc 

where Nmc is the number of second class events in a generated set. 

The main goal of our fits was to find the minimal model compatible with 
our data. To this end, we write the matrix element as follows: 

mpa)\' = + ■ J^r.™ + zx ■ Ji-°°P , 

where X stands for an admixture under study — pa, /ii(1170)7r, p^p^, 
a2(1320)7r or 7r(1300)7r. In the above equation, expressions for Jx""*^ are 
taken from Appendix ^ while complex factors Z are components of the vec- 
tor A. 

2.4 Comparison of 7r+7r~27r° data with simulation 

Since the matrix element for the channel wtt^ has a well-known structure, 
one can use the events from the first class to test the adequacy of the to- 
tal MC simulation ||T^. Figs, show the distributions over Mi„^(7r*7r"), 

^recoil 

(tt^), Mj„^(7r+7r-), M„„(7r°7rO), cos(V'^±^o), cos(-0^+^-), cos{ip^o^o) 
and Mrecoii{T^^) for the events of the first class at the beam energy of 690 
MeV. Points with errors are the data, while the histograms correspond to the 
simulation of the processes wtt^ and ai(1260)7r. Good consistence of the data 
and MC makes us confident that MC simulation adequately reproduces both 
the kinematics of produced particles and the detector response to them. 

Similar distributions for the events of the second class are shown in 
Figs. One can see that the process tt+tt 27r° is satisfactorily described 
in the minimal model in which there are two intermediate states wtt" and 
ai(1260)7r only. Similar consistence is observed at other energies: we have 
also examined the energy points 2Ebeam = (1.28 ± 0.01) and (1.18 ± 0.03) 
GeV. 

To determine the admixture of other possible mechanisms we extend the 
minimal model above by adding each of the other states one by one and 
performing the fit. The results of these fits are shown in Table |^. From 
its third column one can see that the relative fractions rx of the additional 
intermediate state X with respect to the ai(1260)7r are small. This confirms 
our assumption of the ai(1260)7r dominance. 

Because of the theoretical uncertainty of the ai(1260)7r matrix element, 
we do not consider the nonnegligible magnitude of the p~^p~ and 7r(1300)7r 
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Table 2: The results of fit to different models with upper limits 



Model 




rjf, % 


Upper limit, % 


ojtt'^ + aiTT 


1264/891 






WTT^ + aiTT + pa 


1256/891 


2 1+^-2 


4.3 


lutt'^ + aiTT + hlTT 


1263/891 


"•^-0.1 


0.4 


Um'^ + ttiTT + a2TT 


1263/891 




0.8 


UJTT^ + aiTT + Tt'tT 


1250/891 




15. 


ijJTl^ + aiTT + 


1246/891 




7.7 



contributions as significant. Instead, we prefer to set upper limits for the 
fraction of these intermediate states. The values of the likelihood function (^) 
for optimal parameters (Lmm) do not contradict to our expectation based on 
the simulation of wtt" and ai(1260)7r. 

2.5 Comparison of 27r+27r~ data with simulation 

Consider now the 27r+27r^ channel. In this case we don't have W7r° in the 
intermediate state and therefore no additional complicacy to check the as- 
sumption of the ai(1260)7r dominance arises. To this end four-track events 
were selected. All tracks should come from the interaction region: a distance 
from track trajectories to the beam axis should be less than 1 cm and the 
vertex position along the beam axis should be inside ±15 cm. After that 
a kinematic fit was performed assuming the 27r+27r~ hypothesis and events 
with /ndf < 2.5 were selected. Under these conditions 28552 events 
remain in the energy range under study. 

Figure^ shows distributions over Mi„^(7r+7r^), Afi„^(7r='=7r^), Mrecoii{T^^) 
and cos('(/'7r+7r- ) for 27r+27r~case. One can see that the hypothesis of the 
ai(1260)7r dominance does not contradict to the data although one should 
note a slight deviation in the spectrum of the invariant masses of the likesign 
pions. This deviation can be possibly explained by the contribution of the D— 
wave or some final state interaction. It can't be explained by the admixture 
of other possible processes because their fractions compared to ai(1260)7r are 
small (see Table ^ . The study of more complicated cases taking into account 
simultaneously a few intermediate states in addition to W7r° and ai(1260)7r 
is now in progress. 
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Figure 5: Distributions over Mi„i,(7r^7r°), Mrecouij^^)-, Minvij^ 
Minvi'^'^'^'^) for 7r+7r"27r'^ events from the class <M 
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Figure 6: Distributions over cos(V'7r±7rf )j cos{iJjtj+t^-), cosiijj^o^o) and 
MrecoiiiT^'^) for 7r+7r^27r° events from the class (Q) 
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Figure 7: Distributions over Mi„i,(7r^7r°), Mrecouij^^)-, Minvij^ 
Minvi'^'^T^'^) for 7r+7r"27r'^ events from the class (O) 



14 



100 - 




cos(\\f), n n 




30 - 



20 



10 - 



80 



60 



40 



20 



500 




1000 



COS(\\f), Tl\'^ 



Figure 8: Distributions over cos(V'7r±7rf )j (^os{iJjt^+^-), cosiip^o^o) and 
Mrecoiii'^'^) for 7r+7r^27r° events from the class (O) 
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Figure 9: Distributions over Mi„t,(7r+7r ), Mi„t,(7r*7r^), Mrecoui'^^) and 
cos (■(/'7r+ TT - ) for 47r* events 



3 Energy dependence of the cross sections 



To obtain the ai(1260)7r contribution to the total cross section, we subtracted 
from the latter the contribution of the wtt" intermediate state. Strictly speak- 
ing, such a procedure loses its meaning when the interference between wtt" 
and ai(1260)7r is strong. In our case, however, the effects of the interference 
is numerically small (~ 5%) because of the small width of the w meson. 

In order to diminish the influence of background from soft photons, the 
additional cut was applied: cos{ijj) < 0.7 where -0 is the angle between the 
photon direction in the tt*^ rest frame and 7r° momentum. 

Cross sections were calculated from the formulae: 



B{uj^iTi)-{l + 6)-L-e' 



where L is the luminosity, e is the detection efficiency from MC, and (5 is a 
radiative correction calculated according to . 

Nuj and Na^ are the numbers of events due to the corresponding reactions 
{e^e~ —> ujtt'^ and e~^e~ —> oitt respectively) under the assumption that 
there is no interference. The values of and iV^j^ were determined from the 
equations: 

A^i = N^-a + Na,-il-l3), 
N2 = N^-il~a) + Na,-l3, 



where A^i and iV2 are the numbers of events in the first and second classes 
respectively, a is the probability for the wtt" event to enter the first class, and 
P is the probability for the ai(1260)7r event to enter the second class. The 
values of a and /? as a function of energy were taken from MC simulation. 

Figure [l0| shows the cross sections obtained for e^e" wtt" and e+e^ 
7r^7r~27r'' with the W7r° contribution subtracted as explained above. Both 
cross sections rise with energy while the relative fraction of the ai(1260)7r 
increases. 

Figure |ll| shows the total cross section vs energy. Only statistical errors 
are shown. The systematic uncertainty for this channel consists of the con- 
tributions from the uncertainties in the event reconstruction, radiative cor- 
rections and luminosity determination. The overall systematic uncertainty 
was estimated to be ~ 15%. 
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Figure 10: Energy dependence of (T(e+e^ wtt") and (T(e+e^ 7r+7r^27r°) 
where the contribution of wtt" is subtracted 

The cross section measured in our experiment is consistent with the pre- 
vious measurement at OLYA and within systematic errors does not con- 
tradict to the recent result from SND |Q. The value of the cross section 
from all three groups is significantly lower than that from ND [ p3| p^ . Also 
shown above 1.4 GeV are results from Orsay [|l5[ ^ and Frascati p6[ groups. 

To determine the cross section e+e~ 2tt^2tt^ the following criteria 
were added: 

. xVndf < 7 

• Omin > 0.67 

• I Si Pi I < OS- Ebearn 

where Ei, pi are pion energies and momenta before the kinematic recon- 
struction, and Omin is a minimal angle of pions with respect to the beam 
axis. 

Figure |l^ shows the total cross section vs energy. Only statistical er- 
rors are shown. The systematic uncertainty for this channel consists of the 
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Figure 11: Energy dependence of the tt+tt 2Tr^ cross section 
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contributions from the uncertainties in the event reconstruction, radiative 
corrections, selection criteria and luminosity determination. The overall sys- 
tematic uncertainty was estimated to be ~ 15%. 




Figure 12: Energy dependence of the 27r+27r cross section 

The obtained cross section is consistent with the previous measurements 
at OLYA ||2l|, ND (2|, |24| and CMD Q and within systematic errors does 
not contradict to the recent result from SND [ p2|. Also shown above 1.4 GeV 
are results from Orsay ^ |l^ and Frascati ||7|, ^ groups. 

Figure [l^ presents the ratio of the cross sections a{e~^e~ 27r+27r^) and 
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CT(e+e TT+TT 27r") where the contribution of wtt*^ is subtracted. The sohd 
curve shows the theoretical prediction based on the ai(1260)7r dominance. 
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Fi gure 13: Energy dependence of the ratio ^A+g-^^+^-2^o\ ■ The contribu- 
tion of wtt" is subtracted. The sohd curve shows the theoretical prediction 
based on the ai(1260)7r dominance 

In the energy range under consideration the experimentally measured 
ratio is about two. Calculations show that this ratio tends to four at very low 
energies close to the threshold of four pion production when the interference 
effects are maximal. The effects of interference decreases with energy since 
the produced resonances are moving at a high velocity. As a result, the ratio 
falls with energy and tends to unity. 

The importance of the interference between the different diagrams cor- 
responding to permutations of identical pions was first noted in p^ . The 
neglection of interference can lead to wrong conclusions about the contribu- 
tions of different intermediate states to the total cross section. 
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4 Conclusion 



The detailed analysis of the process e+e^ — > 7r+7r~27r'' unambiguously demon- 
strates that the dominant contribution to the cross section comes from lutt^ 
and p^n^TT^ intermediate states whereas the p^Tr^Tr" state is not observed. 
Moreover, the p'^n^n^ state is completely saturated by the ai(1260)7r mech- 
anism. The latter also dominates in the cross section of the process e+e^ 
27r"'"27r~. This conclusion is based on the data sample corresponding to the 
integrated luminosity of 5.8 pb~^ collected with the CMD-2 detector in the 
energy range 1.05 - 1.38 GeV. 

Comparison of the experimental data with the calculations shows that 
account of the interference of different amplitudes with various intermediate 
resonances but the identical final state drastically changes predictions for the 
differential distributions and total cross sections. 

Through the hypothesis of CVC |^ experimentally tested to be valid 
within the 3-5% accuracy , we can relate the values of the four pion cross 
sections measured in e+e^ to the hadronic spectra in corresponding four 
pion decays of the r-lepton. Moreover, our observation of the ai(1260)7r 
dominance should be also applied to r-lepton decays. 

One should note that such a conclusion is in contradiction to that made by 
ARGUS collaboration [|2|. However, their conclusion based on the distribu- 
tion over the dipion mass is not compatible with the recent data from ALEPH 
fl^ . The model used by ARGUS didn't take into account the interference 
of the amplitudes corresponding to various intermediate states whereas our 
observation unambiguously points to significance of the interference effects. 

Our predictions for various two- and three-body invariant mass spectra 
will be addressed in more detail in the forthcoming paper which will com- 
pare them to the existing r-lepton data from ARGUS jl^, CLEO jl^ and 
ALEPH 0. 

The detailed investigation of the energy behavior of exclusive channels 
of e+e^ annihilation into hadrons is crucial for understanding the structure 
of the initial hadronic state created by the virtual photon (/9(1450), p(1700) 
etc.). Our work considers this problem in the energy range below 1.4 GeV. At 
higher energy, such analysis can't be performed because e'^e~ annihilation 
data are not available now. One can hope to move forward by using data on 
semihadronic decays of the r-lepton. 

Our statement about the oi (1260)7r dominance can be used for verification 
of theoretical models in which the problem of intermediate states has been 
discussed. In particular, recently there were numerous speculations about the 
possible existence of vector hybrids claiming that hybrids can be distinguished 
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from the normal qq states by their decays, e.g. in the ai(1260)7r p9| . 
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Appendix 



A Models of substructure in e+e — > 47r 

In this section we discuss the details of the model which has been used for 
description of the experimental data. Since the initial hadron state (referred 
to as p) is created by a virtual photon {e~^e~ — > 7* — > p) and then decays 
into four pions, it has the p-meson quantum numbers {I^J^^ = l"*"! ). 
We use the notation for the polarization vector of this state. 

We assume that the main contribution to the amplitude of the process in 
the energy range under study is given by the intermediate state resonances 
having the masses close to the threshold of pn production (U)-(0)- In the 
case of broad resonances (ai (1260),a2(1320),7r(1300) etc.), the form of their 
propagators 1 /D(q) is very important for analysing the data. We learned that 
it is necessary to take into account the dependence of the imaginary part of 
the propagators (width) on virtuality while the corresponding corrections to 
the real part which can be expressed through the imaginary part by dispersion 
relations are not so important. 

We represent the function D{q) in the form 

Diq)=q^-M^ + zMr^^ (9) 

where M and F are the mass and width of the corresponding particle, and 
the function g{s) describes the dependence of the width on the virtuality. If 
q^ ~ AP then D ~ iMT, in accordance with the usual definition of mass 
and width of the resonance. In the case of the p-meson we used the following 
representation for the function gp{q)'. 

5p(s) = s-V2(,_4m2)3/2. (10) 



A.l The contribution of ai(1260)7r 
The quantum numbers of the ai(1260) resonance are 

ing into account the quantum numbers of the pion, we can write the matrix 
elements corresponding to the processes p{P) — > ai(q)7r(p) and ai{q) — > 
p{P'Mp) as 

T{p -> anr) = Fpa,.e''''\P^,~e, - P,~e^,)q^Al* <j>''* , 
T(ai ^ pn) = Fa.p^e^^'qMPl^t " P'A*W , (H) 
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where a, b, c are isospin indices, A° and ej^ are the polarization vectors of oi- 
and p-mesons, (p^ is the pion wave function, Fa^^p-^ and Fpa^r are form factors 
depending on the virtuahty of initial and final particles. The explicit form 
of these form factors in the energy range considered is not very essential (it 
contributes to the theoretical uncertainty of the model) . The matrix element 
of the transition T{p tttt) reads 

T{p ^ nn) = Fp^^e'^'^elipf , (12) 

where p^f-* and are 4-momenta of the corresponding pions. Due to the 
helicity conservation only transverse space components (with respect to elec- 
tron and positron momenta) of the 4- vector are not zero. We denote them 
as ej^. For unpolarized electrons and positrons a square of matrix element 
for the process p — + 47r is of the form 



T|2=|J^|2. (13) 



Using (11) and ( |12[ ) we obtain the following expression for the contribution 
of the ai(1260)-meson to the current J in the process 
p ^ 7r+(pi)7r+(p2)7i""(p3)7i"^(p4) : 

J^i^ = G[tQi(pi,p2,P3,P4) +tai(pi,P4,P3,J32) +tai(p2,Pl,P3,P4) 
+^at{P2,Pi,P?.,Pl) + tai(pi,P2,P4,P3) + tai (Pl , P3 , P4 , P2) 
+ tai(p2,Pl,P4,P3) +tai(P2,P3,P4,Pl)] , (14) 

where 

tai Pl,P2,P3,P4) = n V7T7 ; V 

DaAP - Pi)Dp[pi +P3) 

•k{{E - £4)[Pi(^^£3 -P4P3) - P3(^^£l -P4Pl)] 

-P4['^l(P4P3) - e3(P4Pl)]} , 

P is the initial 4-momentum { ~ E — 0) ^ pi — (£i,Pj) , 1/D^(g) and 
1/Dp{q) are propagators of the ai- and p- mesons, Fa^q) is the form factor, 
G is the function of the initial energy E proportional to the amplitude of p 
creation. Similarly to (p^), we obtain for the contribution of ai(1260) to the 
current J in the process p 7r"'"(pi)7r~(p4)7r°(p2)7r''(p3): 

J<[r°° = G[tai(pi,p2,P3,P4) -ta,(p4,P2,P3,Pl) (16) 
+ tai(pi,P3,P2,P4) - toi(P4,P3,P2,Pl)] ■ 
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The function (s) in the propagator of ai reads: 

aM = FliQ) j 



£2Pl - eiP2 ^ £2P3 - £3P2 



Dp{j)l+P2) Dp{p2+P3) 

dpi dp2 dps 5 (pi +P2+P3-Q) 



(17) 



2ei2e22£3(27r)5 

where Q° = -/i and Q = 0. As a form factor, we used the function F{q) = 
(l + TO2^/A2)/(l + g2/A2) with A ~ 1 GeV. 

A. 2 The contribution of cun 

The amplitude of the process p{V) — > u{q)Tr{p) has the form: 

T{p ^ wtt) = Fp^^e,,ua0V,,quele% , (18) 

where is the polarization vector of the w-meson. The matrix element of 
the transition lo pn can be written in the similar form. The cj-mcson 
contributes only to a channel p 7r+(pi)7r~(p4)7r''(p2)7''°(p3) . The corre- 
sponding current is equal to 

J+-00 = G^[t^{p2,P4,Pl,P3) -tu;{P2,Pl,P4,P3) 

- twiP2,P3,Pl,P4:)] + {P2 ^ P3) , (19) 

where 

tu,{pi,P2,P3,P4) = n r ^ ^ (20) 

X{('?4P3 - £3P4)(PlP2) - P2(£4PlP3 " £3PlP4) 
-e2[P3(PlP4) - P4(PlP3)]} > 

e, is the energy of the corresponding pion, F^^{q) is the form factor. Since 
the width of the u is small, we set ^^(s) = 1 in the propagator D^{q). 



A. 3 The contribution of /ii(1170)7r 

The quantum numbers of the ft.i(1170) resonance are sim- 
ilarly to the w-meson it gives the contribution to the mixed channel only. The 
matrix elements for the transitions p{P) hi{q)Tr°{p) and hi{q) p{P)Tr{p) 
have the form 

T{p ^ /ii7r°) = Fph^AP^^.^ - P.e^)qpK4,^* , 
T{hi ^ pn) = Fh,p^q^H,{P^el* - P.e';*)^* , (21) 
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where iJ^ is the polarization vector of hi. Using 
the following representation for the /ii-meson contribution to the current J 
in the decay p tt~^ {pi)Tr^ {pi)^'^ {p2)t^^ (ps) ■ 



21) and ([12|) we obtain 



J/1",™ = GhA^hi{Pl,P4.,P3,P2) -th,{p4,Pl,P3,P2) 

- tfti(Pl,P3,P4,P2)] + {P2 ^ Ps) , (22) 



the current t/^ is given by with the change of indices ai — > /ii , and the 
function g^.^ (s) in the propagator of hi is: 

2 



e2Pl - eiP2 e2P3 - e3P2 *^3Pl " ^iPs 



Dp{pi+P2) Dp{p2+P3) Dp{pi+P3) 

dPi dps S'^^'^ {pi +P2+P3-Q) 
2ei2£22£3(27r)5 



(23) 



A. 4 The contribution of p+p 



One more contribution to the mixed channel amplitude comes from the pro- 
cess p p'^ p~ 47r. The matrix element corresponding to the transition 
p{P) p'^{p)p^{q) reads 

T{p^p+p-) = Fpp+p-{P^~e,-P,~e^) (24) 
A{Pt^e+* -Pae+*)(g^e-* - fce^)"*) - (/i ^ v)] 

where and e~ are the polarization vectors of p+ and respectively. 
Using (p4|), we obtain the contribution of p'^ p~ to the current J in the decay 
p ^ 7r+(pi)7r"(p4)7r°(p2)7r°(p3) : 



GppFl{pi+p2)Fl{p3+Pi) 

Dp{pi+p2)Dp{p3+P4) 
x[Pl(e3P2P4 - £4P2P3) - P2(e3PlP4 " e4PiP3) 
-P3(eiP2P4 - e2PiP4) + P4(^lP2P3 " £2PiP3)] + (2^3). 



A.5 The contribution of 7r(1300)7r 

The matrix elements for the transitions p{P) ^ Ti'{q)'ri{p) and 7r'(q) 
p{P)tt{p) are of the form (for brevity tt' = 7r(1300)): 

T{p ^ ttV) = Fp^,^e^''\Ppe, - P,~ep)q^,pA"'* 4>''* , 

r(7r' ^ ptt) = F^,p^e'''''qppAP^el* - F,e^*)0'>^* . (26) 
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Then we get for the contribution of 7r(1300) to the current J in the decay 
p ^ 7r+(pi)7r+(p2)7i""(p3)7i""(p4) : 

J;^'"^ = G^' [t^'bl>-P2,P3,P4) +t^'(p2,Pl,P3,P4) +t^'(pi,P2,P4,P3) 

+t7r'(P2,Pl,P4,P3) + K^{Pi,P2,Pl,Pi) + t^'{p4,P2,Pl,P3) 
+ t^'(?53,Pl,P2,P4) +t^'(p4,Pl,P2,P3)] ■ (27) 

In the mixed channel p 7r+(pi)7r^(p4)7r*^(j52)7i'^(p3) the corresponding cur- 
rent reads: 



-00 



Gtt' [t^'ipi,P2,P3,P4:) + t7r'(Pl'-P3,P2,P4) 

+ t7r'(P4,Pl,P3,P2) +t^'(p4,Pl,P2,P3)] • (28) 



It ( |27| ) and (|2^) we use the notation 

F^'iP-Pi) 



t7r'(Pl,P2,P3,P4) 



D^,{P-pi)Dp{p2 +P4) 

For the function gTr'{s) in the tt' propagator, one has 



ijP2P3-PAP3){p2+P4fVl- (29) 



7.'(s) = F2(g) 



£lP2P3 - £2PlP3 e3PlP2 - i^2PlP3 



Dp{pi +P2) 



Dp{p2 +P3) 



rfPl dp2 rfp3 ^^'^^ (Pl + P2 + P3 ~ Q) 

2£i2e22e3(27r)5 



(30) 



A. 6 The contribution of ap 

The quantum numbers of the <t resonance are 

iGjPC ^ o+0^ 
element of the transition p{P) o'{q)p^{p) is of the form: 

T{p (Tp°) = Fp„p{Ppe^ - P^ep)q^el(l)l. 

The corresponding contribution to the current J in the decay 
p —> 7r+(pi)7r+(p2)7r^(p3)7r~(p4) reads: 

= G„[t„{pi,P2,P3,PA) +t„{p2,Pi,P3,Pi) 
+ t^(pi,P2,P4,P3) + tCT(p2,Pl,P4,P3)] • 



The matrix 



(31) 



(32) 



In the mixed channel p 7r+(pi)7r {pi)'^^(p2)''^'^{P'i) the current has the 
form 

J + -"" = G,W(pi,P2,P3,P4), (33) 
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where 

t^(Pi,P2,P3,P4 = ■ TTTT ■ T £4Pi - eiPi) . 34 

The function ga{s) in the propagator of a is equal to 

= (l-4mVs)'/^ (35) 

A.7 The contribution of a2(1320)7r 

The quantum numbers of the a2(1320) resonance are 

jGjPC = 1-2++. The 

matrix elements for the transitions p{P) — s- a2{q)T^{p) and 02(9) p{P')'K{p) 
can be written in the form 



T{p a27r) = Fpa^T^e " efj,^^pxPfj.et,ppAl*^pjqi^<j)' 

2p7r£ £ fi,vp\PpGjj PpA^^ 



T{a2 pit) = Fa,p^e'''''ep,,pxP'el*PpA\n^cb'* (36) 



where is the polarization tensor of 02 . The contributions of the 02- 
meson to the current J in the processes p 7r+(pi)7r+(p2)7r~(p3)7r~(p4) and 
p 7i'^(pi)7r^(p4)7r°(p2)7r°(p3) are given by the formulae ( p^ ) and ( pj| ) with 
the substitution t^^ — s- ta2 , where 



Da[P - Pi)Dp{pi +P3) 

x{£;(p4 X P2)[pi X P3,P2] + b2P4 - ip2q) {p4.q) / ml^] 

X [Pl(£3P2P4 - '^2P4P3) - P3('^lP4P2) ~ e2P4Pl) 
+P2(eiP3P4 - e3PlP4)]} • 

Here q = P — p^, and P is the initial 4- momentum. 

For the fimction ga2{s) in the propagator of 02, we obtain: 



rfPi rfP2 rfPa 6^'^\pi +P2+ P3 - Q) 

2£i2£22£3(27r)5 



where 

Dp{p2 + P3) P>p{Pl + P3 
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